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Abstract 

The perturbative QCD predictions concerning deep inelastic scattering at low 
X are summarized. The theoretical framework based on the leading log 1/x re- 
summation and kt factorization theorem is described and some recent developments 
concerning the BFKL equation and its generalization are discussed. The QCD ex- 
pectations concerning the small x behaviour of the spin dependent structure func- 
tion gi{x, Q^) are briefly summarized and the importance of the double logarithmic 
terms which sum contributions containing the leading powers of asln?{l/x) is em- 
phasised. The role of studying final states in deep inelastic scattering for revealing 
the details of the underlying dynamics at low x is pointed out and some dedicated 
measurements, like deep inelastic scattering accompanied by an energetic jet, the 
measurement of the transverse energy flow etc., are briefly discussed. 
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Perturbative QCD predicts that several new phenomena will occur when the param- 
eter X specifying the longitudinal momentum fraction of a hadron carried by a parton 
(i.e. by a quark or by a gluon) becomes very small 0] . The main expectation is that 
the gluon densities should strongly grow in this limit, eventually leading to the parton 
saturation effects []T], 0, The small x behaviour of the structure functions is driven 
by the gluon through the g ^ qq transition and the increase of gluon distributions with 
decreasing x implies a similar increase of the deep inelastic lepton - proton scattering 
structure function F2 as the Bjorken parameter x decreases 0. The Bjorken parameter 
X is, as usual, defined as a; = Q^/ {2pq) where p is the proton four momentum, q the four 
momentum transfer between the leptons and = —q^. The recent experimental data are 
consistent with this perturbative QCD prediction that the structure function F2{x,Q'^) 
should strongly grow with the decreasing Bjorken parameter a; |^. 



The growth of structure functions with decreasing parameter x is much stronger than 
that which would follow from the expectations based on the "soft" pomeron exchange 
mechanism with the soft pomeron intercept a so ft ~ 1-08 as determined from the phe- 
nomenological analysis of total hadronic and real photoproduction cross-sections |10[ . 



Small X behaviour of structure functions for fixed reflects the high energy be- 
haviour of the virtual Compton scattering total cross-section with increasing total CM 
energy squared W"^ since W"^ = Q'^{l/x — l). The Regge pole exchange picture would 



therefore appear quite appropriate for the theoretical description of this behaviour. 



The high energy behaviour of the total hadronic and (real) photoproduction cross- 
sections can be economically described by two contributions: an (effective) pomeron with 
its intercept slightly above unity (~ 1.08) and the leading meson Regge trajectories with 
intercept 0^(0) ~ 0.5 |jTO|. The reggeons can be identified as corresponding to p,u!, f 



or A2 exchange (s) depending upon the quantum numbers involved. All these reggeons 
have approximately the same intercept. One refers to the pomeron obtained from the 
phenomenological analysis of hadronic total cross sections as the "soft" pomeron since 
the bulk of the processes building-up the cross sections are low pt (soft) processes. 



The Regge pole model gives the following parametrization of the deep inelastic scat- 
tering structure function F2{x,Q'^) at small x 



F2{x,Q') = j:m' 



}X 



(1) 



The relevant reggeons are those which can couple to two (virtual) photons. The (singlet) 
part of the structure function F2 is controlled at small x by pomeron exchange, while the 
non-singlet part -^2^'^ = F|' — F2 by the A2 reggeon. Neither pomeron nor A2 reggeons 
couple to the spin structure function gi{x,Q'^) which is described at small x by the ex- 
change of reggeons corresponding to axial vector mesons [f 



13 1 i.e. to Ai exchange for 



the non-singlet part gi 



NS 



P 

9i 



gi etc. 



9i 



NS 



(2) 



The reggeons which correspond to axial vector mesons are expected to have very low 
intercept (i.e. < etc.). 
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The high energy behaviour which follows from perturbative QCD is often refered to 
as being related to the "hard" pomeron in contrast to the soft pomeron describing the 
high energy behaviour of hadronic and photoproduction cross-sections. 

In Fig. 1 we summarize the present experimental situation on the variation of the total 
virtual Compton scattering cross-section with total CM energy W for different photon 
virtualities which range from the real photoproduction (Q^ = 0) to the deep inelastic 
region |[T^. The change of high energy behaviour with the scale is evidently present 
in the data. 



The purpose of this talk is to review some of the recent theoretical developments con- 
cerning deep inelastic scattering in the limit of small x. The next section is devoted to 
the discussion of the BFKL equation which resums the leading powers of ln{l/x) and to 
the unified evolution equation |T^, |T^, |T^, |Ty[ based on coherence and angular ordering 
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The kt factorization theorem 21, 22, 



its connection with collinear factorization 



as well as the small x resummation effects within the QCD evolution formalism are also 
discussed in this section. In sec. 3 we present the theoretical expectation for the small x 
behaviour of the spin dependent structure function gi{x,Q'^) concentrating for simplicity 
on its non-singlet part. The novel feature in this case is importance of the double loga- 
rithmic terms i.e. of those terms in the perturbative expansion which contain powers of 
asln^{l/x). Finally in sec. 4 we discuss some of the dedicated measurements like deep 
inelastic + jet events and transverse energy flow etc. which aim at revealing the dynam- 
ical details of small x physics. Sec. 5 contains a brief summary and outlook. 



2. The BFKL pomeron and QCD predictions for the 
small X behaviour of the structure function F2 

At small X the dominant role is played by the gluons and the basic dynamical quantity is 
the unintegrated gluon distribution f{x, Qf) where x denotes the momentum fraction of 
a parent hadron carried by a gluon and Qt its transverse momentum. The unintegrated 
distribution /(a;, Qf) is related in the following way to the more familiar scale dependent 
gluon distribution g{x, Q"^): 



xg{x,Q^ 



Q 



(3) 



In the leading ln{l/x) approximation the unintegrated distribution f{x,Qf) satisfies the 
BFKL equation ^, |l6l which has the following form: 



1 da^ 
x' 



{q + Q 



/(x,Q?) = r(x,Q?)+ 

-f(x',{q + Q,r)-fix',Q',)eiQ', 



where 



a. 



3a= 



vr 



(4) 



(5) 



The first and the second terms on the right hand side of eq. correspond to real gluon 
emission with q being the transverse momentum of the emitted gluon, and to the virtual 
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corrections respectively. p{x,Ql) is a suitably defined inhomogeneous term. 

After resumming the virtual corrections and " unresolvable" gluon emissions (g^ < 
jj?) where /i is the resolution defining the "resolvable" radiation, equation @) can be 
rearranged into the following "folded" form: 

f{x,Ql) = f%x,Ql)+ 
where A/j which screens the 1/z singularity is given by: 



a. 



2 2;' J^t2 



and 

/«,.,Q?,^/;f^A„(£0?)f|g^ (8) 



c//n(l/x') 

Equation (|^) sums the ladder diagrams (see Fig. 2) with the reggeized gluon exchange 
along the chain with the gluon trajectory adQt) = 1 — ^[ln{Ql / ^i^) . 

For the fixed coupling case eq. (^) can be solved analytically and the leading behaviour 
of its solution at small x is given by the following expression: 

fix, Ql) ~ {Ql) 2 exp - : y ; 9 

with 

Xbfkl = 4/n(2)a, (10) 

A" = «,28C(3) (11) 

where the Riemann zeta function ({3) ~ 1.202. The parameter Q is of nonperturbative 
origin. 



The quantity 1 + Xbfkl is equal to the intercept of the so - called BFKL pomeron. Its 
potentially large magnitude (~ 1.5) should be contrasted with the intercept agoft ~ 1-08 
of the (effective) "soft" pomeron which has been determined from the phenomenological 
analysis of the high energy behaviour of hadronic and photoproduction total cross-sections 



The solution of the BFKL equation reflects its diffusion pattern which is the direct 
consequence of the absence of transverse momentum ordering along the gluon chain. 
In this respect the BFKL dynamics is different from that based on the (leading order) 
Altarelli-Parisi evolution which corresponds to strongly - ordered transverse momenta. 
The interrelation between the diffusion of transverse momenta towards both the infrared 
and ultraviolet regions and the increase of gluon distributions with decreasing x is a 
characteristic property of QCD at low x. It has important consequences for the structure 
of the hadronic final state in deep inelastic scattering at small x. 
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In practice one introduces the running coupling as{Q^) in the BFKL equation (^). 
This requires the introduction of an infrared cut-off to prevent entering the infrared region 
where the couphng becomes large. The effective intercept Xbfkl found by numerically 
solving the equation depends weakly on the magnitude of this cut-off p5| . The impact of 
the momentum cut-offs on the solution of the BFKL equation has also been discussed in 



refs. [E6|, 1271 



The solution of the BFKL equation is obtained most directly by solving the cor- 
responding equation for the moment function 



dx 



10 X 
^2\ _ rOf 



a 



2 



[q + Qt 



f{u;,Qi) = f{uj,Qi) + 



(12) 



(13) 



This equation can be diagonalised by a Mellin transform. The solution for the Mellin 
transform /(a;, 7) of the moment function f{uj,Q^) is: 



1 - t^(7) 



(14) 



where 



ir(7) = 2^(1) - ^(7) - ^(1 -7) (15) 

is the Mellin transform of the kernel of eq. (p^ . The function "^{z) is the logarithmic 
derivative of the Euler F function. The Mellin transform f{uj,'y) is defined by: 



/(^,7) 



dQ 
Ql 



(16) 



and hence the function f{uj, Qf) is related to /(cj, 7) through the inverse Mellin transform 



1 /•1/2+ioo 

TT- 1 dj{Q',yf{u,^). 

ZTTI Jl/2-ioo 



(17) 



The poles of f{uj,'y) in the 7 plane define the anomalous dimensions of the moment 
function f{uj,Qf) ||2^. The (leading twist) anomalous dimension 'jgg{u,as) of f{uj,Q'f) 
gives the following behaviour of f{u,Qf) at large 



where 



R{as,uj) 



^IggK^iCH-s) \-y=-ygg{uj,as) 



-I -1 



(19) 



The anomalous dimension will also, of course, control the large behaviour of the 
moment function g{u, Q^) of the integrated gluon distribution 



Qt 



(20) 
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which has the following form: 

-)2\ laai'^AB) 



-g{uj,Q^)=R{as,u)f{uj)[^J (21) 



where we have introduced the moment function of the input distribution 

fiu) = Pico, J = 7,,(u;,«,))(g^)7..(-,a.). (22) 

Equations (0)and (|2^) follow directly from equations (|18]), ([T9|) and (|20|) . It may be 
seen from eq. ( pT] ) that the BFKL singularity affects the " starting" gluon distribution at 
Q2 = g2 through the factor R |2|. 

It follows from eq.(|l5D that the anomalous dimension 'jggiuj, as)) is the solution of the 
following equation: 

-^K{jgg{u;,as)) = l. (23) 
uj 

It is a function of only one variable — i.e. Jggiuj, cxs) lgg{~)- The solution of eq. ( p3D 
makes it possible to obtain the anomalous dimension 7c,c,(^) as a power series of ^ 0] 

7..(% = i:^n(|)". (24) 

UJ „=i V / 

This power series corresponds to the leading ln{l/z) expansion of the splitting function 



■Pgg\^y ^s) 



^ . X ^ (a,/n(l/z))"-i 
zPggiz, as) = asY. , ^ (25) 



which controls the evolution of the gluon distribution. 



The structure functions F2^l{x,Q'^) are described at small x by the diagram of Fig.3 
which gives the following relation between the structure functions and the unintegrated 
distribution /: 

F2A^, Q') = f'^f ^F^Li^'. Ql Q')f{-, QD- (26) 

The functions F2°£{x', Ql-,Q'^) may be regarded as the structure functions of the off-shell 
gluons with virtuality Qf. They are described by the quark box (and crossed box) dia- 
gram contributions to the photon-gluon interaction in the upper part of the diagram of 
Fig. 3. The small x behaviour of the structure functions reflects the small z {z = x/x') 
behaviour of the gluon distribution f{z,Qf). 



Equation (|26|) is an example of the " kt factorization theorem" which relates measurable 
quantities (like DIS structure functions) to the convolution in both longitudinal as well as 
in transverse momenta of the universal gluon distribution f{z, Qf) with the cross-section 
(or structure function) describing the interaction of the "off-shell" gluon with the hard 



probe p2|, ^ . The kt factorization theorem is the basic tool for calculating the observable 
quantities in the small x region in terms of the (unintegrated) gluon distribution / which 
is the solution of the BFKL equation. 
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The leading - twist part of the kf factorization formula can be rewritten in a coUinear 
factorization form. The leading small x effects are then automatically resummed in the 
splitting functions and in the coefficient functions. The kt factorization theorem can in 
fact be used as the tool for calculating these quantities. Thus, for instance, the moment 
function Pgg{uj,as) of the splitting function is represented in the following form (in the 
DIS scheme): 

p , , 7|,(!f)^r(c^ = 0,7 = 7..(!f)) 

Pgg{U},as) = 9 (27) 

where F2°^{uj,'y) is the Mellin transform of the moment function F2°^ {u , Qf , Q"^) i.e. 



, ,1/2+ioo ^ 

Ql Q') = —:l .^ . d^F'2°%u, 7) ( ^ (2^ 



ZTTZ Jl/2-ioo \Qt 

Representation (^7]) generates the following expansion of the splitting function Pgg{z, as) 
at small z: 

zPgg{z, a,) = ^zP^'\z) + {asf J2 (29) 

The first term on the right hand side of eq. (|29|) vanishes at z = 0. It should be noted that 
the splitting function Pgg is formally non-leading at small z when compared with the split- 
ting function Pgg . For moderately small values of z however, when the first few terms in 
the expansions ( |2^ ) and ( P^ dominate, the BFKL effects can be much more important in 
Pgg than in Pgg. This comes from the fact that in the expansion (^) all coefficients 6„ are 
different from zero while in eq. ( ^^ we have C2 = C3 = The small x resummation 



effects within the conventional QCD evolution formalism have recently been discussed in 
refs. [p8| , p9l , pO| , pT| . In Fig. 4 we show the results of this analysis for the structure 



function F2{x,Q'^) |2^ where the starting parton distributions have been assumed to be 
non-singular. In this figure we compare the standard two-loop results (dashed curves) 
with those which contain the leading ln{l/x) resummations. The dot-dashed curve cor- 
responds to keeping the small x corrections only in the gluon anomalous dimension while 
the solid and dotted curves show the effect of including the small x resummation in the 
splitting function Pgg as well. (The two curves were obtained by applying two different 
prescriptions of implementing the momentum sum rule which is violated in the LLl/x 
approximation). One can see from this figure that at the moderately small values of x 
which are relevant for the HERA measurements, the small x resummation effects in the 
splitting function Pgg have a much stronger impact on F2 than the small x resummation 
in the splitting function Pgg. This reflects the fact, which has already been mentioned 
above, that in the expansion (|29| ) all coefficients 6„ are different from zero while in eq. 
we have 02 = 03 = 0. 



A more general treatment of the gluon ladder than that which follows from the BFKL 
formalism is provided by the CCFM equation based on angular ordering along the gluon 



chain |T^, |T8|. This equation embodies both the BFKL equation at small x and the 
conventional Altarelli-Parisi evolution at large x. The unintegrated gluon distribution / 
now acquires dependence upon an additional scale Q which specifies the maximal angle of 
gluon emission. The CCFM equation has a form analogous to that of the "folded" BFKL 
equation (||): 

f{x,QlQ') = f%x,QlQ')+ 
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q + Qt 



(30) 



where the theta function Q{Q — qx/x') reflects the angular ordering constraint on the 
emitted gluon. The "non-Sudakov" form-factor Aij(2;, , q^) is now given by the following 
formula: 



^R{z,Qt,q'^) = exp 



J z z J q 



(31) 



Eg . (1301) still contains only the singular term of the g ^ gg splitting function at small z. 
Its generalization which would include remaining parts of this vertex (as well as quarks) 



is possible. The numerical analysis of this equation was presented in ref. |[18 



In Fig. 5 we show the results for the structure function F2 calculated from the kt 



factorization theorem using the function / obtained from the CCFM equation p2[. We 
confront these predictions with the most recent data from the HI and ZEUS collabora- 
tions at HERA 1^, as well as with the results of the analysis which was based on the 
Altarelli-Parisi equation alone without the small x resummation effects being included in 



the formalism 0. In the latter case the singular small x behaviour of the gluon and 
sea quark distributions has to be introduced in the parametrization of the starting distri- 
butions at the moderately large reference scale = Ql (i.e. Ql ~ AGeV"^ or so) ||3^. One 



can also generate steep behaviour dynamically starting from non-singular "valence-like" 

)2 = 0.35Ge\/2 n 



parton distributions at some very low scale Ql 
gluon and sea quark distributions exhibit "double logarithmic behaviour" 



In the latter case the 



F2(x, Q') ~ exp 2Ji{Q\ Ql)ln{l/x) 



where 



i{Q\Ql) 



dq^ 3^5 (g^ 



(32) 



(33) 



For very small values of the scale Ql the evolution length ^{Q"^, Qq) can become large for 
moderate and large values of and the "double logarithmic" behaviour ( p2D is, within 
the limited region of x, similar to that corresponding to the power like increase of the 
type x~^, A ~ 0.3. This explains similarity between the theoretical curves presented in 
Fig. 5. The theoretical results also show that an inclusive quantity like F2 is not the best 
discriminator for revealing the dynamical details at low x. One may however hope that 
this can be provided by studying the structure of the hadronic final state in deep inelastic 
scattering and this possibility will be briefly discussed in sec. 4. 

Although the BFKL and CCFM equations for the gluon distribution become equiva- 
lent in the leading ln{l/x) approximation they begin to differ if less inclusive quantities 
are considered. In the case of the observables which are not infrared safe (like multiplicity 
of jets etc.) the CCFM formalism generates double logarithmic ln{l/x) terms i.e. terms 
which contain powers of asln'^{l/x). This is of course closely related to the fact that 
the angular ordering constraint serves also as the infrared regulator and so the transverse 
momentum integrations generate additional powers of ln{l/x) besides the "conventional" 
ones which come from the integration over the longitudinal momenta. Also the virtual 
corrections contain double logarithmic terms as well. This can easily be seen from the 



definition of the "non-Sudakov" form factor (plf). The double logarithmic terms exactly 
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cancel between the real emission contributions and virtual corrections to the gluon dis- 
tribution but remain present in the less inclusive quantities. On the other hand in the 
BFKL case only the single logarithmic terms (i.e. powers of asln{l/x)) are present in 
all quantities but the observables which are not infrared safe diverge for /z^ — where fi 
is the resolution parameter. The gluon distribution itself remains of course finite in this 
limit. 



Several new interesting results have been obtained in the formal studies of the high 
energy (or small x) limit in QCD which go beyond the leading logarithmic approximation 
0, ^ ^ . The important theoretical tool in this case is the effective field theory 

where the basic objects are the reggeized gluons and the effective action of this effective 
field theory obeys conformal invariance [BBl E^. Theoretical analysis simplifies in the 



large Nc limit. One can discuss both the pomeron which appears as the bound state of 
two (reggeized) gluons, and the odderon (i.e. the bound state of three reggeized gluons), 
as well as bound states of many reggeized gluons [^]. Conformal invariance is also very 



useful for analysing the BFKL pomeron away from the forward direction 40 1 as well 
as the the triple pomeron and more complicated vertices 



The genuine next-to-leading ln(l/a;) corrections to the BFKL equation can be present 
in all relevant quantities i.e. in the particle-particle-reggeon vertex, the reggeon-reggeon- 
particle vertex and in the gluon Regge trajectory (The reggeon here corresponds 

to the reggeized gluon). Besides that one has also to include additional region of phase- 
space which goes beyond strong ordering of longitudinal momenta. 



It should finally be emphasised that in impact parameter representation the BFKL 
equation offers an interesting interpretation in terms of colour dipoles P2|. 



4. Small X behaviour of the nonsignlet unpolarized 
and polarized structure functions 

The discussion presented in the previous Section concerned the small x behaviour of the 
singlet structure function which was driven by the gluon through the g ^ qq transition. 
The gluons of course decouple from the non-singlet channel and the mechanism of gener- 
ating the small x behaviour in this case is different. 

The simple Regge pole exchange model predicts in this case that 

^(x, Q') = Fiix, Q') - F,^ix, Q') ~ x'-^-^^'^ (34) 

where 0^2(0) is the intercept of the A2 Regge trajectory. For 0:^2(0) ~ 1/2 this behaviour 
is stable against leading order QCD evolution. This follows from the fact that the leading 
singularity of the moment 'Jqq{uj) of the splitting function Pqq{z): 

iH = I -z-P,,{z) (35) 

is located at u; = and so the (nonperturbative) A2 Regge pole at a; = 0^2(0) ~ 1/2 
remains the leading singularity controlling the small x behaviour of the non-singlet struc- 
ture function. 
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The novel feature of the non-singlet channel is the appearence of the double log- 
arithmic terms i.e. powers of asln^il/x) at each order of the perturbative expansion 
p5| , |i6| , 0, ^ . These double logarithmic terms are generated by the ladder diagrams 
with quark (antiquark) exchange along the chain. The ladder diagrams can acquire cor- 
rections from the " bremsstrahlung" contributions H^, EU| which do not vanish for the 
polarized structure function g^^{x,Q'^) P9| . 



In the approximation where the leading double logarithmic terms are generated by 
ladder diagrams illustrated in Fig. 6 the unintegrated non-singlet quark distribution 
f^^{x, k1) [q^^ = u + u — d — d) satisfies the following integral equation : 



cNS 



1 dz dQ 



z JQl Q 



/2 

12 Jq 



X 



1/2 \ 



where 



-a. 



(36) 



(37) 



and Ql is the infrared cut-off parameter. The unintegrated distribution f^^{x, Qf) is, as 
usual, related in the following way to the scale dependent (nonsinglet) quark distribution 



^^^{x,Q'y. 



Q 



(3^ 



The upper limit Q'^/z in the integral equation ( |36D follows from the requirement that 
the virtuality of the quark at the end of the chain is dominated hj Q^. A possible non- 
perturbative A2 reggeon contribution has to be introduced in the driving term i.e. 



f7(a:,Q,^)~a; 



-aAo (0) 



(39) 



at small x. 



Equation (|36|) implies the following equation for the moment function f^^{uj,Ql 



a. 



Equation ( ^OD follows from (^) after taking into account the following relation: 

r ^^Kl - ^) ^ M^^^^ - ^^^^ ^ (i)^ ^^^^^ - 



(40) 



(41) 



For fixed coupling otg equation (^OD can be solved analytically. Assuming for simplicity 
that the inhomogeneous term is independent of Ql (i.e. that f^^{LJ, Ql) = C{uj) ) we get 
the following solution of eq. (^OD : 



f'^'iuj,Q-t) = C{uj)R{a,,uj) 



2\7 ias,t^) 



where 



7 l«s,^J 



(42) 



(43) 
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and 

R{as,uj) = z ■ (44) 

Equation (|43|) defines the anomalous dimension of the moment of the non-singlet quark 
distribution in which the double logarithmic ln{l/x) terms i.e. the powers of ^ have 
been resummed to all orders. It can be seen from ( ^31) that this anomalous dimension has 
a (square root) branch point singularity aX u = u where 

uj = 2v^. (45) 

This singularity will of course be also present in the moment function f^^{uj, Ql) itself. It 
should be noted that in contrast to the BFKL singularity whose position above unity was 
proportional to a^, is proportional to ^/a's - this being the straightforward consequence 
of the fact that equation (^0]) sums double logarithmic terms (%)". This singularity gives 
the following contribution to the non-singlet quark distribution f^^{x,Q^) at small x: 

For small values of the QCD coupling this contribution remains non - leading in compar- 
ison to the contribution of the A2 Regge pole. 

As has been mentioned above the corresponding integral equation which resums the 
double logarithmic terms in the spin dependent quark distributions is more complicated 
than the simple ladder equation (^) due to non-vanishing contributions coming from 
bremsstrahlung diagrams. It may however be shown that , at least as far as the non- 
singlet structure function is concerned, these contributions give only a relatively small 
correction to u. The main interest in applying the QCD evolution equations to study the 
spin structure function is that the naive Regge pole expectations based on the exchange of 
low-lying Regge trajectories become unstable against the QCD perturbative "corrections". 
The relevant reggeon which contributes to gi^{x, Q^) is the Ai exchange which is expected 
to have a very low intercept aAi(O) < 0. The perturbative singularity generated by the 
double logarithmic ln{l/x) resummation can therefore become much more important 
than in the case of the unpolarized case where it is hidden behind leading A2 exchange 



contribution. Even if we restrict ourselves to leading order QCD evolution ^ then 
the non-singular x~°'^i^^^ behaviour (with q;^^(0) < ) becomes unstable as well and the 
polarized quark densities acquire singular behaviour: 



Ag^^(x, Q2) ^ exp{2^^^S{Q^)ln{l/x)) (47) 



where 



This follows from the fact that APqq{z) = Pqq{z) where Pqq{z) and APqg(z) are the split- 
ting functions describing the evolution of the spin independent and spin dependent quark 
distributions respectively and from the fact that Pqq{z) —>■ const as z ^ 0. 



e'iQ') = J %<is{q') (48) 



The introduction of the running coupling effects in equation (^OD turns the branch 
point singularity into the series of poles which accumulate at a; = |H6| . The numerical 
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analysis of the corresponding integral equation, with the running coupling effects taken 
into account, gives an effective slope , 



dln{l/x) 



(49) 



with magnitude A(x, Q^) ~ 0.2 — 0.3 at small x [^]. The result of this estimate sug- 
gest that a reasonable extrapolation of the (non-singlet) polarized quark densities would 
be to assume an x~^ behaviour with A ~ 0.2 — 0.3. Similar extrapolations of the spin- 
dependent quark distributions towards the small x region have been assumed in several 
recent parametrizations of parton densities |53|, |54|, pSl ^ 



The small x behaviour of the 
spin dependent structure function gi has also been disscussed in refs. EB . 



5. The structure of the hadronic final state in deep 
inelastic scattering at low x 

It is expected that absence of transverse momentum ordering along the gluon chain, which 
leads to the correlation between the increase of the structure function with decreasing x 
and the diffusion of transverse momentum should reflect itself in the behaviour of less 
inclusive quantities than the structure function F2{x,Q'^). The dedicated measurements 
of low X physics which are particularly sensitive to this correlation are the deep inelas- 
tic plus jet events, transverse energy flow in deep inelastic scattering, production of jets 
separated by the large rapidity gap and dijet production in deep inelastic scattering. The 
diagrammatic illustration of these measurements is presented in Fig. 7. 

In principle deep inelastic lepton scattering containing a measured jet can provide 



a very clear test of the BFKL dynamics at low x I^O], ^ |6j. The idea is to study 
deep inelastic (x, Q^) events which contain an identified jet (xj, k'^j^^) where x « Xj and 
~ k^j. Since we choose events with Q"^ f» k^j the leading order QCD evolution (from 
k^j to Q^) is neutralized and attention is focussed on the small x, or rather small x/xj 
behaviour. The small x/xj behaviour of jet production is generated by the gluon radiation 
as shown in the diagram of Fig. 7a. Choosing the configuration fa k^j we eliminate by 
definition gluon emission which corresponds to strongly ordered transverse momenta i.e. 
that emission which is responsible for the LO QCD evolution. The measurement of jet 
production in this configuration may therefore test more directly the (x/xj)"'^ behaviour 
which is generated by the BFKL equation where the transverse momenta are not ordered. 
The recent HI results concerning deep inelastic plus jest events are consistent with the 
increase of the cross-section with decreasing x as predicted by the BFKL dynamics . 

A conceptually similar process is that of the two-jet production, with the jets separated 
by a large rapidity gap Ay, in hadronic collisions or in photoproduction as illustrated in 
Fig. 7c Besides the characteristic exp{XAy) dependence of the two-jet cross- 



section one expects significant weakening of the azimuthal back-to-back correlations of 
the two jets. This is the direct consequence of the absence of transverse momentum or- 
dering along the gluon chain in the diagram of Fig. 7c. 

Another measurement which should be sensitive to the QCD dynamics at small x is 
that of the transverse energy flow in deep inelastic lepton scattering in the central re- 
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gion away from the current jet and from the proton remnant as illustrated in Fig. 
p. 



7b. 



The BFKL dynamics predicts in this case a substantial amount of transverse en- 
ergy which should increase with decreasing x. The experimental data are consistent with 
this theoretical expectation |Q. Absence of transverse momentum ordering also implies 
weakening of the back-to-back azimuthal correlation of dijets produced close to the pho- 
ton fragmentation region (see Fig. 7d) 



3c 



6. Summary and conclusions 

In this talk we have briefly described the QCD expectations for deep inelastic lepton 
scattering at low x which follow from the BFKL dynamics. It leads to the indefinite 
increase of gluon distributions with decreasing x which is correlated with the diffusion of 
transverse momenta. This increase of gluon distribution implies a similar increase of the 
structure functions through g qq transitions. Besides discussing the theoretical and 
phenomenological issues related to the description of the structure function F2 at low x 
we have also emphasised the role of studying the hadronic final state in deep inelastic 
scattering. 

The indefinite growth of parton distributions cannot go on forever and has to be even- 
tually stopped by parton screening which leads to the parton saturation. Most probably, 
however, the saturation limit is still irrelevant for the small x region which is now being 
probed at HERA. 

We have also discussed the small x behaviour of the spin dependent structure func- 
tion Qi focussing for simplicity on its non-singlet component which, at small x , can be 
(approximately) described by ladder diagrams with the quark (antiquark) exchange. The 
novel feature in this case is the appearance of double logarithmic terms i.e. contributions 
which contain powers of agln^il/x). The perturbative QCD effects become significantly 
amplified for the singlet spin structure function due to mixing with the gluons. The simple 
ladder equation may not however be applicable for an accurate description of the double 
logarithmic terms in the polarized gluon distribution AG. 

We have limited ourselves to the large region where perturbative QCD is expected 
to be applicable. Specific problems of the low Q^, low x region are discussed in ref. ||70| . 
Finally let us point out that the change of the dynamics with the relevant scale is clearly 
visible in the data (see Fig. 1) and its satisfactory explanation is perhaps one of the most 
challenging problems of to-day. 



Acknowledgments 

I thank Robin Devenish and Mike Whalley for organizing an excellent meeting. I thank 
Barbara Badelek, Krzysztof Golec-Biernat, Alan Martin and Peter Sutton for most en- 
joyable research collaborations on the problems presented in this lecture. I am grateful 
to Grey College and to Physics Department of the University of Durham for their warm 
hospitality. This research has been supported in part by the Polish State Committee for 



13 



Scientific Research grant 2 P03B 231 08 and the EU under contracts nO. CHRX-CT92- 
0004/CT93-357. 



References 

[1] L.N. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1. 

[2] B. Badelek et aL, Rev. Mod. Phys. 64 (1992) 927. 

[3] A.D. Martin, Acta. Phys. Polon, B25 (1994) 265. 

[4] J. Kwiecihski, Nucl. Phys. B (Proc. Suppl.) 39 B,C (1995) 58. 

[5] E.A. Kuraev, L.N.Lipatov and V.S. Fadin, Zh. Eksp. Teor. Fiz. 72 (1977) 373 
(Sov. Phys. JETP 45 (1977) 199); Ya. Ya. Bahtzkij and L.N. Lipatov, Yad. Fiz. 
28 (1978) 1597 (Sov. J. Nucl. Phys. 28 (1978) 822); J.B. Bronzan and R.L. Sugar, 
Phys. Rev. D17 (1978) 585; T. Jaroszewicz, Acta. Phys. Polon. Ell (1980) 965. 

[6] L.N. Lipatov, in "Pcrturbative QCD", edited by A.H. Mueller, (World Scientific, 
Singapore, 1989), p. 441. 

[7] A.J. Askew et al., Phys. Rev. D47 (1993) 3775; Phys. Rev. D49 (1994) 4402. 

[8] HI collaboration: A. de Roeck et al., preliminary measurements to be pubhshed 
in the Proc. of the Workshop on DIS and QCD, Paris, 1995, DESY preprint 95 - 
152. 

[9] ZEUS collaboration: B. Foster, to be published in the Proc. of the Workshop on 
DIS and QCD, Paris, 1995, DESY preprint 95-193. 

[10] A. Donnachie and P.V. Landshoff, Phys. Lett. B296 (1992) 257. 

[11] P.D.B. Collins, "An Introduction to Rcggc Theory and High Energy Physics", 
Cambridge University Press, Cambridge, 1977. 

[12] B.L. loffe, V.A. Khoze and L.N. Lipatov, "Hard Processes", North Holland, 
Amsterdam-Oxford- New York- Tokyo , 1984. 

[13] J. EUis and M. Karliner, Phys. Lett. B231 (1988) 72. 

[14] A. Donnachie and P.V. Landshoff, Z. Phys. C61 (1994) 161. 

[15] A. Levy, these proceedings. 

[16] M. Ciafaloni, Nucl. Phys. B296 (1988) 49. 

[17] S. Catani, F. Fiorani and G. Marchesini, Phys. Lett. B234 (1990) 339; Nucl. 
Phys. B336 (1990) 18; G. Marchesini, in Proceedings of the Workshop "QCD at 
200 TeV", Erice, Italy, 1990, edited by. L. CifareUi and Yu. L. Dokshitzer (Plenum 
Press, New York, 1992), p. 183; G. Marchesini, Nucl. Phys. B445 (1995) 49. 

[18] J. Kwiecihski, A.D. Martin, P.J. Sutton, Phys. Rev. D52 (1995) 1445. 



14 



[19] B. Andersson, G. Gustafson and J. Samuelsson, Lund preprint LU TP 95-13. 
[20] Yu.L. Dokshitzer et al., Rev. Mod. Phys. 60 (1988) 373. 

[21] S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B242 (1990) 97; Nucl. 
Phys. B366 (1991) 657; J.C. Colhns and R.K. Elhs, NucL Phys. B 360 (1991) 3; 
S. Catani and F. Hautmann, Nucl. Phys. B427 (1994) 475. 

[22] M. Ciafaloni, Phys. Lett. 356 (1995) 74. 

[23] J. Kwiecihski, A.D. Martin, Phys. Lett. B353 (1995) 123. 

[24] T. Jaroszewicz, Phys. Lett. B116 (1982) 291. 

[25] A.D. Martin, J. Kwiecihski and P.J. Sutton, Nucl. Phys. B (Proc. Suppl.) A29 
(1992) 67. 

[26] J.C. Collins, P.V. Landshoff, Phys. Lett. B276 (1992) 196. 

[27] J.F. McDermott, J.R. Forshaw and G.G. Ross, Phys. Lett. B349 (1995) 189. 

[28] R.K. Ellis, Z. Kunszt and E.M. Levin, Nucl. Phys. B420 (1994) 517; Erratum-ibid. 
B433 (1995) 498. 

[29] R.K. Ellis, F. Hautmann and B.R. Webber, Phys. Lett. B348 (1995) 582. 
[30] R.D. BaU and S. Forte, Phys. Lett. 351 (1995) 313. 

[31] J.R. Forshaw, R.G. Roberts and R.S. Thorne, Phys. Lett. B356 (1995) 79. 

[32] J. Kwiecihski, A.D.Martin and P.J. Sutton, Durham preprint DTP/95/94. 

[33] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D50 (1994) 6734; Phys. 
Lett. 354 (1995) 155. 

[34] M. Gliick, E. Reya and A. Vogt, Z.Phys. C67 (1995) 433. 

[35] A. de Rujula et al., Phys. Rev. DIO (1974) 1649. 

[36] L.N. Lipatov, Nucl. Phys. B452 (1995) 969; talk given at the International Con- 
ference on Elastic and Diffractive Scattering, "Frontiers in Strong Interactions", 
Chateau de Blois, France, June 1995. 

[37] R. Kirschner, L.N. Lipatov and L. Szymanowski, Nucl. Phys. B425 (1994) 579; 
Phys. Rev. D51 (1995) 838. 

[38] L.D. Fadeev and G.P. Korchemsky, Phys. Lett. B342 (1995) 311; G.P. Korchem- 
sky, to be pubhshed in the Proc. of the Workshop on DIS and QCD; G.P. Korchem- 
sky, Stony Brook preprint ITP-SB-95-25; R. Janik, Cracow preprint TP JU- 18-95. 

[39] A.H. Mueller and W.K. Tang, Phys. Lett. B284 (1992). 

[40] J. Bartels et al., Phys. Lett. B348 (1995) 589. 

[41] J. Bartels, M. Wiisthoff, Z. Phys. C66 (1995) 157; J.. Bartels, L.N. Lipatov and 
M. Wiisthoff, DESY preprint 95-171. 



15 



[42] A.H. Mueller, Nucl. Phys., B415 (1994) 373; A.H. Mueller and B. Patel, Nucl. 
Phys. B425 (1994) 471; A.H. Mueller, Nucl. Phys. B437 (1995) 107; Chen Zhang 
and A.H. Mueller, Nucl. Phys. B451 (1995) 579. 

[43] V.S. Fadin, talk given at the International Conference on Elastic and DifFractive 
Scattering, "Frontiers in Strong Interactions", Chateau de Blois, France, June 
1995. 

[44] A.R. White, talk given at the International Conference on Elastic and Diffractive 
Scattering, "Frontiers in Strong Interactions", Chateau de Blois, France, June 
1995. 

[45] V.G. Gorshkov et al, Yad. Fiz. 6 (1967) 129 (Sov. J. Nucl. Phys. 6 (1968) 95); 
L.N. Lipatov, Zh. Eksp. Teor. Fiz. 54 (1968) 1520 (Sov. Phys. JETP 27 (1968) 
814). 

[46] J. Kwiecihski, Phys. Rev. D26 (1982) 3293. 

[47] R. Kirschner and L.N. Lipatov, Nucl. Phys. B213 (1983) 122. 

[48] B.I. Ermolaev, S.I. Manayenkov and M.G. Ryskin, DESY preprint 95-017. 

[49] J. Bartels, B.I. Ermolaev and M.G. Ryskin, DESY preprint 95 - 124. 

[50] M.A.Ahmed and G.G. Ross, Phys. Lett. B56 (1975) 385; Nucl. Phys. Bll (1976) 
298; G. Altarelh and G. Parisi, Nucl. Phys. B126 (1977) 298. 

[51] R.D. Ball, S. Forte and G. Ridolfi, Nucl. Phys. B444 (1995) 287. 

[52] J. Kwiecihski, Durham preprint DTP 95/98 (to appear in Acta Physica Polonica). 

[53] C. Bourrely and J. Soffer, Phys. Rev. D51 (1995) 2108; Nucl. Phys. B445 (1995) 
341; Marseille preprint CPT-95/P.3224. 

[54] M. Gliick, E. Reya and W. Vogelsang, Phys. Lett. B359 (1995) 201. 

[55] T. Gehrmann and W.J. Stirhng, Durham preprint DTP/95/62. 

[56] J. Bliimlein and A. Vogt, DESY preprint 95-175. 

[57] S.D. Bass, P.V. Landshoff, Phys. Lett. B336 (1994) 537. 

[58] F.E. Close, R.G. Roberts, Phys. Lett. B336 (1994) 257. 

[59] A.D. Martin, review talk presented at the Blois Workshop "The Heart of the 
Matter", June 1994, Blois, France. 

[60] A.H. Mueller, J. Phys. G17 (1991) 1443. 

[61] J. Bartels, M. Loewe and A. DeRoeck, Z. Phys. C54 (1992) 635. 
[62] W.K.Tang, Phys. Lett. B278 (1992) 363. 

[63] J. Kwiecihski, A.D. Martin and P.J. Sutton, Phys. Rev. D46 (1992) 921; Phys. 
Lett. B287 (1992) 254. 



16 



[64] 



HI Collaboration, DESY preprint 95 - 108. 



[65] V. del Duca, Phys. Rev. D49 (1994) 4510. 
[66] W.J. Stirling, Nucl. Phys. B423 (1994) 56. 

[67] J. Kwiecinski, A.D. Martin, P.J. Sutton and K.Golec-Biernat, Phys. Rev. D50 
(1994) 217; K.Golec-Biernat, J. Kwiecinski, A.D. Martin and P.J. Sutton, Phys. 
Lett. B335 (1994) 220. 

[68] J.R. Forshaw and R.G. Roberts, RAL preprint 94-0228. 

[69] A.J. Askew et al., Phys. Lett. B338 (1994) 92. 

[70] B. Badelek and J. Kwiecinski, Warsaw preprint IFD/1/94 (to appear in Rev. Mod. 
Phys.). 

Figure captions 

1. The total virtual photon cross-section plotted as a function of for different values 
of the photon virtuality Q^. The curves correspond to the theoretical parametriza- 
tions [0, [1^. In particular the skewed curve is the GRV prediction for = 
0.356*6^^. The figure is taken from ref. []15 |. 



2. Diagrammatic representation of the BFKL equation (^. 

3. Diagrammatic representation of the kt factorization formula (0). 

4. Theoretical predictions for the structure function F2{x,Q'^) based on the Altarelli- 
Parisi evolution equations with the leading ln{l/x) terms resummed in both the Pgg 
and Pgg splitting functions (solid and dotted curves where the dotted curve is based 
on the different prescription for imposing the momentum sum rule). The dot-dashed 
curves show predictions when the leading ln{l/x) resummation is included only in 
the splitting function Pgg. The dashed curves correspond to the two- loop prediction. 
The theoretical curves are confronted with the 1993 data from HERA. The figure 
is taken from ref. [|2D|. 



A comparison of the HERA measurements of F2 [§1^ with the predictions based 
on the kt factorization formula ( ^Bf ) using for the unintegrated gluon distributions / 
the solutions of the CCFM equation (^) (continuous curve) and of the approximate 
form of this equation corresponding to setting Q{Q — q) in place of Q{Q — qx/x') 
and = 1 (dotted curve). We also show the values of F2 obtained from coUinear 



factorization using MRS (A') [|33| and GRV ||3^ partons (the figure is taken from ref. 



H). 



6. The ladder diagram with quark (antiquark) exchange. 

7. Diagrammatic representation of the processes testing BFKL dynamics, (a) Deep 
inelastic scattering with the forward jet. (b) flow in deep inelastic scattering, 
(c) Production of jets separated by the large rapidity gap Ay. (d) Dijet production 
in deep inelastic scattering (the figure is taken from ref. []59| ). 
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